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High-precision electroweak experiments such as parity-violating M0ller scattering can provide 
indirect access to physics at multi-TeV scales and play an important complementary role to the 
LHC research program. However, before physics of interest can be extracted from experimental 
data, electroweak radiative corrections, which can significantly reduce the cross-section asymme- 
try, must be calculated with an unprecedented completeness and accuracy. Although the two-loop 
corrections are strongly suppressed relative to the one-loop corrections, they can no longer be dis- 
missed for the upcoming precision experiments. We evaluate a full gauge-invariant set of one-loop 
and several types of two-loop radiative corrections for the parity-violating e~e~ —¥ e~e~ (y)(yy) 
scattering asymmetry by combining two distinct but mutually-reinforcing techniques: semiauto- 
matic, precise, with FeynArts and FormCalc as base languages, and by hand, with some approx- 
imations. For 1 1 GeV relevant for the ultra-precise MOLLER experiment planned at JLab, the 
results obtained by two approaches are in excellent agreement, which gives us assurance that our 
calculations are error-free. 
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1. Introduction 

Polarized M0ller scattering has been a well-studied low-energy reaction for close to eight 
decades now []l|], but recently has attracted especially active interest from both experimental and 
theoretical communities as an excellent tool for measuring parity- violating (PV) asymmetries 

The first observation of parity violation in M0ller scattering was made by the E-158 ex- 
periment at SLAC which studied M0ller scattering of 45- to 48-GeV polarized electrons 
on the unpolarized electrons in a hydrogen target. Its result at low Q 2 = 0.026 GeV 2 , Apy = 
(1.31 ±0.14 (stat.)±0.10 (syst.)) x 10~ 7 |@] allowed one of the most important parameters in 
the Standard Model (SM) - the sine of the Weinberg angle - to be determined with an accuracy 
of 0.5% (sin dw = 0.2403 ± 0.0013 in the MS scheme). A very promising experiment measur- 
ing the ep scattering asymmetry soon to be completed at JLab, Qweak aims to determine 
sin 2 dw with relative precision of 0.3%. The next-generation experiment to study e~e~-scattering, 
MOLLER (Measurement Of a Lepton Lepton Electroweak Reaction), planned at JLab follow- 
ing the 11 GeV upgrade, will offer a new level of sensitivity and measure the PV asymmetry 
in the scattering of longitudinally polarized electrons off an unpolarized target to a precision of 
0.73 ppb. That would allow a determination of the weak mixing angle with an uncertainty of 
±0.00026 (stat.) ±0.00013 (syst.) [|0, or about 0.1%, an improvement of a factor of five in frac- 
tional precision when compared with the E-158 measurement. At such precision, any inconsistency 
with the Standard Model will signal new physics, so M0ller scattering experiments can provide in- 
direct access to physics at multi-TeV scales and play an important complementary role to the LHC 
research program 

Obviously, before we can extract reliable information from the experimental data, it is neces- 
sary to take into account higher-order effects, i.e. electroweak radiative corrections (EWC). The 
inclusion of EWC is an indispensable part of any modern experiment, but will be of the paramount 
importance for MOLLER. Of course, a significant theoretical effort has been dedicated to one-loop 
radiative corrections already. In [||], we gave a short review of the literature and calculated a full 
gauge-invariant set of the one-loop EWC both numerically, with no simplifications, and by hand, 
in an approximate but compact form. The total correction was found to be close to —70%, and we 
found no significant theoretical uncertainty coming from the largest possible source, the hadronic 
contributions to the vacuum polarization. The dependence on other uncertain input parameters, like 
the mass of the Higgs boson, was below 0.1%. 

A much larger theoretical uncertainty in the prediction for the asymmetry may come from the 
two-loop corrections. According to 0, the higher-order corrections are suppressed by a factor of 
up to 5%, depending on a type of corrections, relative to the one-loop result. However, since the 
one-loop weak corrections for the M0ller scattering are so large and since the 11 GeV MOLLER 
experiment is striving for such an unprecedented precision, theoretical predictions for its scattering 
asymmetry must include not only a full treatment of one-loop radiative corrections (NLO) but also 
leading two-loop correction (NNLO). 

One way to find some indication of a size of the higher-order contributions is to compare 
results obtained with different renormalization schemes. Our calculations [|K]] in the on-shell and 
CDR (Constrained Differential Renormalization) schemes show a difference of about 11%, which 



is comparable with the difference of 10% between MS [11] and the on-shell scheme [Bp. It is also 
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worth noting that although two-loop corrections to the cross section may seem to be small, it is 
much harder to estimate their scale and behaviour for such a complicated observable as the PV 
asymmetry to be measured by the MOLLER experiment. 

We approach a formidable task of calcutating the NNLO EWC in stages. The two-loop 
EWC to the Born (~ MqMq) cross section can be divided into two classes: the <2-part induced 
by quadratic one-loop amplitudes (~ M\M^), and the 7 -part corresponding to the interference of 
the Born and two-loop diagrams (~ 2ReM 2 M^). The goal of this paper is to calculate the <2-part 
and a sub-set of the T-part including the gauge invariant set of boson self energies and vertices of 
two-loop amplitude M 2 . Although it is too early to draw any final conclusions, we believe that the 
two-loop EWC may be larger that previously thought. 



2. General Notations and Matrix Elements 

The cross section of polarized M0ller scattering with the Born kinematics: 

e- (£1) + e~ -> e~ (k 2 ) + e~ (p 2 ), (2.1) 

can be expressed as: 

3 3 

a = —\M + M { + M 2 \ 2 » —(M M+ + 2ReM l M+ +M { M+ + 2ReM 2 M+), (2.2) 
2s 2s 

where a = do /d cos 6 and 6 is the scattering angle of the detected electron with 4-momentum k 2 
in the center-of-mass system of the initial electrons. The 4-momenta of initial (ki and p\) and final 
(k 2 and p 2 ) electrons generate a standard set of Mandelstam variables: 

s = {k l + Pl f, t = (k l -k 2 )\ u = {k 2 - Pl f. (2.3) 

Mo, Mi and M 2 are the Born (ff(a)), one-loop {{?(a 2 )) and two-loop (^(a 3 )) amplitudes (matrix 
elements), respectively. 

The one-loop amplitude M\ as a sum of boson self-energy (BSE), vertex (Ver) and box di- 
agrams. We use the on-shell renormalization scheme from [12, 13], so there are no contribu- 
tions from the electron self-energies. The question of the dependence of EWC on renormalization 
schemes and renormalization conditions (within the same scheme) was addressed in our earlier 
paper [10]. Now we present the one-loop complex amplitude as the sum of IR and IR-finite parts 
Mi = +m{. The IR-finite part m{ can be found in [Q and for the IR part we have: 

M\ = %- 8? Mo , 8f = 45 log 4= , (2-4) 

27C y/S 

where X is the photon mass and the complex value B can be presented in the following form (see, 
for example, [|15|]): 

B = \og-^--l-i%. (2.5) 

m z s 

i f 

Analogously, the two-loop amplitude is the sum M 2 = M^ + M 2 , where 

M| = ^5M+^(f) 2 W 2 M . (2.6) 
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Note that the structure of first term in (2.6) is the same as in (2.4) in terms of the soft photon 
factorization. 

Now we should make sure that the infrared divergences are cancelled. In a similar way used 
for the amplitudes, we present the differential cross sections as sums of A -dependent (IRD-terms) 
and A -independent (infrared-finite) parts: Oi = + a{ , Oq t = Oq t + Oq t . The one-loop cross 
section is already carefully evaluated with full control of the uncertainties in [§]. The simplest form 
for IRD-terms are: 



a 



af = ^a , o* = (£) W+2Re(W) o\ = Re MY + 28? 8 



7T 1 u ' ^ \2n 



r _A 



/ a \ 2 



V27T/ 



a . 



(2.7) 

The imaginary part of the total cross section cancels out in the sum Q- and r -parts due to follow- 
ing properties: 5*5** + Re (8?) 2 = 2(Re 8?) 2 , and Re(5 1 / 5 1 A *) + Re(8{S?) = Re(8{) Re(8?). 
Thus, in the following sections we can ignore the imaginary part, i. e. 8? — > Re5j and B — > ReB. 

3. Bremsstrahlung and Cancellation of Infrared Divergences 

To evaluate the cross section induced by the emission of one soft photon with energy less then 
(0, we follow the methods of Jl6| ] (see also [0]). Then, this cross section can be expressed as: 
o r = oj + cj, where oJ 2 have the similar factorized structure based on the factorization of the 
soft-photon bremsstrahlung: o\ 2 = ^ [— 8? + Ri] Gb,i, where 

Rl = _4Slog ^ - (log - 2 - l) 2 + 1 - t- + log 2 -. (3.1) 
2(0 \ m z / 3 f 

The first part of the soft-photon cross section, oj, cancels the IRD at the one-loop order, while the 
second part, o 2 , cancels the IRD at the two-loop order, with half of o 2 going to the cancellation of 
the IRD in the <2-part and the other half going to treat IRD in the T-part. At last, the cross section 



induced by the emission of two soft photons with a total energy less then co is calculated in [ 14 ] as 



ayy= l -{^)\{-8^+R l f-R 2 )a Q , (3.2) 



where \ is a statistical factor and R 2 = §7T 2 5 2 . 

Combining all the terms, we get the infrared-finite result at both the first and second orders: 

a , 



a NL0 = ffi + o\ = - [Rx + 8{}o°, (3.3) 



K 



tfNNLO = o^ + Oj +a 2 r + a rr ': 



^) 2 [ Rl 8{ + l -R]- l -R 2 + 8^ + 8(}a° 



= o J + o J B + o J Q + o^ (3.4) 

where 



^ = ^iOnlo, oj 1 = - l -(^) 2 (R 2 l +R 2 )o i \ (3.5) 
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4. Numerical Results and Conclusions 



For the numerical calculations at the central kinematic point of MOLLER (£^=11 GeV, 



6 = n/2) we use a, mw, niz and lepton masses as input parameters in accordance with [18]. The ef- 



fective quark masses which we use for the vector boson self-energy loop contributions are extracted 
from shifts in the fine structure constant due to hadronic vacuum polarization AaPjl(m§)=0. 02757 



J19Q. For the mass of the Higgs boson, we take mh = 125 GeV and for the maximum soft photon 
energy we use a> = 0.05-^, according to [[| and [|2C|]. 

Let us define the relative corrections to the Born cross section due to a specific type of contri- 
butions (labeled by C) as 

8 C = (o c -o°)/o°, C = NLO,0,B,2,r,NNLO. 



In the text below the term 'T-part" corresponds to the contributions of a gauge invariant set of the 
BSE and vertices only. The parity-violating asymmetiy is defined in a traditional way: 



Oll + Olr - o RL - o RR 



<5lL - Orr 



OLL + C LR + Or L + Orr LL + 20 LR + Orr ' 

and the relative corrections to the Born asymmetry due to C-contribution are defined as 



(4.1) 



(Alr ~A lr )/A lr . 



Our numerical estimations are presented in the table below, with dots in 4-th and 7-th columns 
denoting the result written in previous line: 



c 


8 C 




C 


8 C 




C 


S c 




NLO 


-0.1144 


-0.6932 


NLO 


-0.1144 


-0.6932 


NLO 


-0.1144 


-0.6932 





0.0457 


0.2347 


...+0 


-0.0687 


-0.3956 


...+(0+B)/2 


-0.1339 


-0.5671 


B 


-0.0848 





...+B 


-0.1535 


-0.4353 


-+Q 


-0.1150 


-0.6392 


Q 


0.0189 


-0.0731 


-+Q 


-0.1345 


-0.5118 


...+(0 + B)/2 


-0.1345 


-0.5118 


T 


0.0119 


-0.1063 


...+T 


-0.1226 


-0.6274 


...+T 


-0.1226 


-0.6274 



As one can see from our numerical data, at the MOLLER kinematic conditions, the part of the 
NNLO EWC we considered in this work can increase the asymmetry by up to ~ 7%. The Q- and T- 
parts do not cancel each other but, on the contrary, are adding up to increase the physical PV effect. 
Clearly, the large size of the investigated parts demands a detailed and consistent consideration 
of the rest of the T-part, which will be the next task of our group. Since the problem of EWC 
for the M0ller scattering asymmetiy is rather involved, a tuned step-by-step comparison between 
different calculation approaches is essential. To make sure that our calculations are error-free, 
we control our results by comparing the data obtained from the equations derived by hand with the 
numerical data obtained with a semi-automatic approach based on FeynArts, FormCalc, LoopTools 
and Form. These base languages have already been successfully employed in similar projects (jSj, 



[ 10]), so we are highly confident in their reliability. In the future, we plan to address the remaining 
two-loop electroweak corrections up to a level required by the plannedl precision of the MOLLER 
experiment and the possible future experiments at ILC. 
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